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We study a spinless two-band model at half-filling in the limit of infinite dimen- 
sions. The ground state of this model in the non-interacting limit is a band- insulator. 
We identify transitions to a metal and to a charge-Mott insulator, using a combina- 
tion of analytical. Quantum Monte Carlo, and zero temperature recursion methods. 
The metallic phase is a non-Fermi liquid state with algebraic local correlation func- 
tions with universal exponents over a range of parameters. 
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Strong correlations can lead to low energy behavior qualitatively different from that 
predicted by band theory. A classic example is provided by the Mott transition where 
a system, which is expected to be a band metal, turns into an insulator as correlations 
become greater than a critical value. Recently, the complementary question of whether 
interactions can drive a band insulator into a metal has been raised in connection with the 
possible (semi)-metallic nature of the "Kondo insulators". From a theoretical point of 
view, it is likely that a metal which exists only as the result of interactions can not be 
analytically continued to a free electron metallic system and hence should not be described 
by Fermi liquid theory. The related question of the possible existence of non-Fermi liquid 
metallic states has been of great recent interest, due to the anomalous normal state properties 
of the high-Tc oxides [|,|. 

In this letter, we study a spinless model at half-filling which displays a band insulating 
ground state in the absence of interactions. We developed two numerical algorithms to 
study the the effects of interactions in infinite dimensions, supplementing the numerical 
results with analytical arguments. We find that, in spite of the simplicity of the model, 
the interactions give rise to a novel correlation-induced insulator (to which we refer as a 
charge- Mott insulator), and a non- Fermi liquid metallic phase with scale invariant low- 
energy correlation functions. In the latter, we identify two parameter regimes, in which the 
associated exponents are continuously varying and universal, respectively. 

We study the following Hamiltonian 

H = ^(ek - /i)cj,Ck + J2(^d - fJ')4di 

k i 

+Y.t{dU + h.c.) + - \Md., - i), (1) 

i i 

where and S are creation operators for two species of spinless electrons. The c-electrons 
have a dispersion e^, while the d-band is dispersionless with local energy e^. These two 
species are coupled through the hybridization, t, and the density-density interaction, V . 
Below we focus on the half-filled case, with = ec = /i = 0. 

We solve this model in the limit of infinite dimensions |^ in which the hopping matrix 
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element tij must be rescaled with \j\fd. The one-particle Green's function, G(k, = 
-/dre^"^"^ < T,V'c.(k,r)V'i(k,r) >, with ^'^ = {d\ c^), can be calculated from the Dyson's 
equation 



G ^(k,iujn) 



iuJn + fJ'-ea t 

t iuJn + A* — ek 



(2) 



The self-energy is momentum independent and can thus be calculated from the local 

Green function, < ip'^ip >Simp' associated with the impurity action 0, 



Simp (G^o) 



dTdT'ij^{T)G~\r - r')ij{r') + f drVic"^ c d --) (3) 
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according to the Dyson equation, S = G~^— < ip'^ip >Simp5 where Go obeys the self- 
consistency equation J2]iG{k,iu!n)/Nsite =< >Simp- already imphcit in the choice 
of the local action @, we will focus on states with no long range (CDW) order 0. 

The effective-impurity action (|^) can be equivalently written in terms of a resonant-level 
model. 



Hi^p = H„ + t{cU + Sc) + (e° - ii)dU + V{nd - -)( 



(4) 



where Ho describes the local c-electron coupled to a non-interacting electron bath with a 
hopping matrix determined from the self-consistency equation. For simplicity we study the 
case of a semi-circular conduction-electron density of states, Po(e) = Z]k^(^ ~ ^i^l^site = 



[2/7rD)y 1 — (e/Z))2, for which the self-consistency equation reduces to the form: 

{G-^)cc = iuJr, + fx- &^G,, (5) 



We developed two numerical techniques to analyze Eqs. 



. The first algorithm is 



a generalization of the method used in previous studies of the Hubbard model. [pHlO|| An 
initial guess for Go is used in the impurity action Eq. @, which is solved by Monte Carlo 
method. Our algorithm is based on a modification of that due to Hirsch and Fye JIT 



and uses a different set of Hubbard- Stratonovich variables [12| and a single determinant. 
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The calculated Green's function is then fed into the self-consistency condition and the 
procedure is iterated until convergence is achieved. 

The second algorithm is based on the recursion method. |T3| , p!4| The one particle Green's 
function can be written as G{z) = G^{z) + G^{z) where G<{z) = J dep±{e)/{z — e) and 
p±{e) = for e^O. We use the continued-fraction expression for G^{z), 

G>{z) = , (6) 

z — al 



r ^ 

z — a2 



and a similar one for G^{z), with a[ and 6[ replaced by a[ and b[. This representation obeys 
the correct analyticity requirements {bl+b{ = 1). It allows us to visualize the self-consistency 
equation, {G~^)cc = i^n — (cc ^ /^) ~ (-D/2)^G'cc5 the inverse of the local Green's function 
< c^c > of the non-interacting Hamiltonian, 



N N-l 



Ho = (ec - p)c^c+^^Jh{{c^cir + h.c.) + ^ a^cl^Q,. + ^ y^6^(ct,Ci+ir + h.c.) 



TV N-\ 



+ + h.c) + a\c\cii + ^h\^-^{c\ci+ii + h.c.) (7) 

^ i=l i=l 

where N is the level of the continued-fraction expression which is infinite in the exact 
representation but taken as finite in practice. The problem can now be reformulated in 
terms of the self-consistent determination of the variables a^, a-, and h\ ||15|. We start 



with an initial guess for these parameters, and solve the ground state |0 > of the resonant- 
level Hamiltonian given by (^) and (0) using the standard sparse matrix diagonalization 
methods. |]14| The new parameters a[, a-, and h\ for the next iteration is then obtained 



from the matrix elements of an orthogonal basis which tridiagonalizes Himp. Specifically, 

< =< /;iHi^pi/; > / < f:\fi >, and hi =< f:\fi > I < f:_,\f:_, >, bi =< >, 

where |/[ >= c|0 >, \f^ >= Hi,^p|/[ > -a^|/[ > and for n > 2, >= Hi^p|/; > 

—al\f^ > —b1\f^_i >. Similar equations, with \f[ >= c^O >, give rise to a[ and b\. 
Compared to the finite temperature Monte Carlo method, this algorithm works directly at 
zero temperature, and is less expansive numerically. 



Using both methods, we find that the model displays two kinds of insulating phases 
separated by a metallic phase. The phase diagram in the V — t parameter space is shown 
schematically in Fig. |1]. To illustrate the phase diagram, we plot in Fig. ^ the d-electron 
local Green's funct ion, Gddi 

function of the Matsubara frequency for two parameters 
with a same value of the hybridization, t = OAD, and different values of the interaction, 
V = —O.ID and V = -~2.5D, respectively. They clearly display the gap in their respective 
single particle excitation spectrum, and also demonstrate the excellent agreement between 
the zero temperature results from the recursion method and the low (but finite) temperature 
results from the Quantum Monte Carlo. Shown in Fig. |^ are Gdd, as well as the c-electron 
local Green's function, Gcc, for the hybridization identical to that used in Fig. 0, t = OAD, 
and an interaction strength in between those used in Fig. ^, = —0.9D. The single-particle 
excitation spectrum is gapless in this case. 

In order to clarify the nature of the phases and the structure of the phase diagram, we 
now consider various limiting cases where analytical arguments can be made. In the non- 
interacting limit (V = 0) S = and the Green function (0) displays two branches of poles, 
-^k — (^k ± \/ (ck)^ + 4t^)/2, describing the fully occupied bonding band (— ) and the empty 
anti-bonding band (+). Thus, a hybridization gap Aq = \ADM-4t^ — D, develops, leading 
to band-insulating behavior. This state becomes unstable whenever the effect of interactions 
is strong enough to close the hybridization gap. 

This can be most easily seen in the atomic limit, D = 0, in which case the behavior of the 
effective impurity model (@-§|) is determined by the local empty, doubly occupied, bonding 
and anti-bonding singly occupied electron states, respectively, 

\E > 
\D > 
\B > 

\A > 

The corresponding eigenvalues are e{E) = e{D) = V"/4, e{B) = —V/i — t, and e{A) = 
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|iV = >= 10 > 



|iV = 2 >= Sc^O > 
1 



\N = 1,6 >= 



V2 



(dt-ct)|0> 



\N = l,ab >= + c^)\0 >, 

v2 



(8) 



—V/A + 1. Three parameter regimes occur in this hmit: (i) For —V < 2t the local bonding 
state, \B >, has the lowest energy, and is separated by a finite gap from lowest single particle 
excitation, leading to band-insulating behavior; (ii) At —V = 2t three states, \B >, \D >, 
and \E >, become degenerate and metallic behavior with gapless single-particle excitations 
sets in; (iii) Finally, for —V > 2t the low-energy charge doublet, \D > and \E >, is separated 
by a gap from the lowest single-particle excitation leading to a situation analogous to the 
Mott insulator of the atomic limit of the half-filled repulsive Hubbard model, in which case 
the spin doublet is the lowest energy state. ^ We thus refer to the resulting phase as a 



charge- Mott insulator. [It 



The stability of these three phases and the position of the boundaries between them, 
can also be understood analytically in the limit of small hybridization. In that case, the 
instability of the band insulator, which occurs through the closing of the hybridization 
gap, is most easily seen in terms of the irrelevance of the hybridization coupling (in the 
sense of the renormalization group) for a sufficiently large value of the attractive inter- 
action, V. The value of the critical coupling, Vd, follows from a scaling analysis of the 
resonant-level model (^, in the vicinity of the phase boundary from the metallic {t ~ 0) 
side. It is obtained from the condition, ej = 1, on the the anomalous dimension of the 
hybridization operator, tt, defined from the long time behavior of the auto-correlation 
function, x(r) = — < T^S {t)c{t)c' {0)d{0) >~ r~^^*. Such a condition implies fl^ 



p{Vci)Vci = -{2/n)tan[Ti{y/2-l)/2]. Here, p{V) = {2/nD)/:Jl - {V/DY is the bath density 
of states at the Fermi level which is itself dependent on \^ as a result of the self-consistency 
requirement (|]). [0] This argument leads to Vd = —0.6D. Upon further increasing the 
attraction beyond —Vd the charge-Mott gap develops at a second critical value, —Vc2- The 



latter can be calculated from the exact solution of the Falicov-Kimball model (at d = oo) p9 
corresponding to the Hamiltonian (|I]) with t = 0, and is identified as the value of V at which 
the local conduction electron density of states at the Fermi energy, Pc(0), starts to vanish. 



For the semi-circular density of states, Pc(0) = {2/txD)^1 — (^)^, ||2^ leading to Vc2 = —D. 
Our numerical results are consistent with the analytical estimates for Vd and V^2- 
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To further characterize the nature of the metalhc phase, we have analyzed numerically 
the low frequency behavior of the local Green's functions. Consider first the limit of small hy- 
bridizations. In this regime, we find that Gd ~ ~i\uJn\~'^ sgnujn- This self-similar behavior is 
consistent with the irrelevant hybridization , which implies a zero renormalized Fermi energy. 
| T7| In this case, the exponent is expected to to be given by the orthogonality exponent, {^Y-, 



with the phase shift determined from the standard formula, 5 = {2/TT)tan^^{7Tp(y)V/2). The 
numerically determined exponents are consistent with this analytical expression. Finally, the 
calculated c-electron local Green's function, Gcc, is nonsingular, as expected from a picture 
of the metallic state in which the local degrees of freedom asymptotically decouple from the 
conduction-electron bath at low energies. 

The most interesting regime of the metallic phase occurs for intermediate to large values 
of the hybridization (as well as the interaction strength). In this regime, our numerical 
results show that, Gdd continues to have a self-similar behavior, with an exponent a close 
to 1, independent of the interaction strength. At the same time, Gcc is again nonsingular, 
suggesting that, even in the intermediate to large hybridization limit (for finite bandwidth, 
2D), a partial decoupling of the d and c electron degrees of freedom occurs at low energies. 

Such a partial decoupling of the local degrees of freedom, and the interaction- 
independence of the exponents, can be understood as a result of the self-consistent modi- 
fication of the electron bath for the effective impurity model using the atomic limit metal- 
lic state as a starting point. In that case, projecting onto the three degenerate low- 
energy atomic states, \B >, \D >, and \E >, in (||), leads to the representation of local 
fermions, = —^{Xbe + ^db) and = -^{^be — ^db) in terms of projection oper- 
ators, Xai3 = \a >< (3\. The corresponding low-energy spectral functions take the form, 
Pc{uj) = Pd{uj) = ^6{uj). In the presence of the hopping term D, these atomic configura- 
tions are coupled to the self-consistent electron bath through the Ho term in (^) which we 
parametrize as. Ho = Z]fc^fc(c^% + H.c.) + J2k^kVkVk- The parameters, ik and Vk, are de- 
termined from the self-consistency condition, (l/NsUe) J2k ^k^i^ ~ ^fe) = (-D/2)^Pc(e), which 
corresponds to taking the imaginary part of Eq. (^. The form of Pc in the atomic limit 



implies that the electron bath 77^ corresponds to a single level, created by r]l, at zero energy. 
Hence, the impurity problem reduces to Hjmp = ~^^[i-^BE + Xdb)Vo + h.c], which acts 
on a Hilbert space spanned by \E >, \B >,\D > and ril\E >, ril\B >,ril\D >. It is easily 
seen that the ground state manifold of this Hamiltonian is doubly degenerate and consists 
of the states = 1 >= (l/v^)(|5 > +ril\E >), and |iV = 2 >= (l/v^)(|D > +ril\B >). 
Recalculating the Green functions of local fermions we obtain a zero energy pole in Gdd and 
pairs of non-zero energy poles, symmetrically centered around zero, in Gcc- This form of 
Gcc then implies that the electron bath for the next iteration is composed of pairs of levels 
at energies ±ej. It can be shown, following the same procedure, that for each pair of levels 
at energies ±e coupled to the local c-electron, there again exist degenerate ground states, 
leading to a zero energy pole in Gad and pairs of non-zero energy poles in Gcc- As the 
iteration proceeds, the poles in Gcc spread over a range of the order of the bare bandwidth 
while a zero-energy pole persists in Gdd, leading to the self-consistent solution with a regular 
Gcc and a self-similar Gdd with the exponent a = 1. 

In summary, we reported a) novel transitions from a band-insulator to a metal and 
to a charge-Mott insulator within a simple spinless two-band model; and b) a novel non- 
Fermi liquid metallic regime characterized by an algebraic local correlation functions with 
universal exponents. These local correlation functions are shown to be associated with a 
non-trivial effective impurity problem. And the very existence of this new non-Fermi liquid 
is closely related to a level-crossing in the atomic limit of the associated impurity problem, a 
benchmark for a class of non-Fermi liquid states in a variety of strongly correlated impurity 



models p^j23|JT7|j2^ . Related ideas have been proposed in a phenomenological discussion 
of strongly correlated lattice models in which non-Fermi liquid behavior arises due to the 
coexistence of single particle states and dispersionless three-body resonances at the chemical 



potential. The results of this letter are expected to have important implications for 
models describing realistic strongly correlated systems such as the "Kondo-Insulators" and 
high Tc copper oxides. In this regard we note that, while the metallic phase in our model 
occurs in a range of attractive interactions, these interactions should be understood as 
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effective parameters derived from the bare interactions in a more complex set of local orbitals 



in real materials. [^] We also note that, in the spinful cases and away from half- filling, there 
can exist an intermediate non-Fermi liquid metallic state with coherent spin excitations and 



incoherent charge excitations. Our results here should describe the charge sector of the 
intermediate phase. These and other related questions are currently under investigation 
within more realistic models. 
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FIGURES 

FIG. 1. The qualitative phase diagram. Here "BI", "M", and "CMI" correspond to the band 
insulating, metallic, and charge-Mott insulating phases respectively. 

FIG. 2. Gdd vs. Matsubara frequency ujn for V = —O.ID and t = OAD (band insulator) 
calculated from the recursion method (solid line) using 10 sites to truncate the impurity Hamilto- 
nian and 100 levels for the continued-fraction expansion, and from quantum Monte Carlo (circles) 
with P = 32. Similar results are shown (dashed line and squares) for V = —2.5D and t = OAD 
(charge-Mott insulator) . 

FIG. 3. Gdd vs. Matsubara frequency Un for V = —0.9D and t = OAD (metal) calculated from 
the recursion method (solid line) with 10 sites and 100 levels and quantum Monte Carlo (circles) 
with (3 = 64. Inset plots the similar results for Gcc- 
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